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ABSTRACT 

We investigate the gravitational instability of galactic discs, treating stars and cold in- 
terstellar gas as two distinct components, and taking into account the phenomenology 
of turbulence in the interstellar medium (ISM), i.e. the Larson-type scaling relations 
observed in the molecular and atomic gas. Besides deriving general properties of such 
systems, we analyse a large sample of galaxies from The Hi Nearby Galaxy Survey 
(THINGS), and show in detail how interstellar turbulence affects disc instability in 
star-forming spirals. We find that turbulence has a significant effect on both the inner 
and the outer regions of the disc. In particular, it drives the inner gas disc to a regime 
of transition between two instability phases and makes the outer disc more prone to 
star-dominated instabilities. 

Key words: instabilities - turbulence - ISM: kinematics and dynamics - ISM: struc- 
ture - galaxies: ISM - galaxies: kinematics and dynamics. 



1 INTRODUCTION 

Today, 30 years after the pioneering work by Larson (1981), 
observations and simulations of the interstellar medium 
(ISM) are revealing its turbulent nature with higher and 
higher fidelity (see, e.g., Elmegreen & Scalo 2004; McKee & 
Ostriker 2007; Romeo et al. 2010). A fundamental aspect of 
ISM turbulence is the existence of scaling relations between 
the mass column density (E), the ID velocity dispersion 
(a), and the size of the region over which such quantities 
are measured (£): 



E oc £ a , a oc £ 



(1) 



The values of a and b, and the range of scales spanned by £ 
depend on which ISM component we consider. In this paper 
we focus on cold interstellar gas, which is highly supersonic 
and hence strongly compressible, and which is known to play 
an important role in the gravitational instability of galactic 
discs (e.g., Lin & Shu 1966; Jog & Solomon 1984a, b; Bertin 
& Romeo 1988, and references therein). 

In the molecular gas, H2, the scaling exponents are 
a w and b « ~, and Eq. (1) holds up to scales of a few 100 
pc. In fact, both Galactic and extragalactic giant molecular 
clouds (GMCs) are fairly well described by Larson's scaling 
laws, E = constant and a oc £ 1 ' 2 , although the uncertain- 
ties are still large (e.g., Larson 1981; Solomon et al. 1987; 
Bolatto et al. 2008; Heyer et al. 2009; Hughes et al. 2010; 
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KaufTmann et al. 2010; Lombardi et al. 2010; Sanchez et al. 
2010; Azimlu & Fich 2011; Ballesteros-Paredes et al. 2011; 
Field et al. 2011; Kritsuk & Norman 2011; Roman-Duval et 
al. 2011; Beaumont et al. 2012). Besides, Larson-type scaling 
relations have now been observed, for the first time, in the 
dense star-forming clumps of a high-redshift galaxy (Swin- 
bank et al. 2011). 

In the atomic gas, Hi, the scaling exponents are instead 
a ~ i and 6 ~ |, and Eq. (1) seems to hold up to scales of 
a few kpc. A Kolmogorov scaling for both <j and E is sug- 
gested by the observed power spectra of H I intensity fluctu- 
ations, and is also consistent with other measurements (e.g., 
Lazarian & Pogosyan 2000; Elmegreen et al. 2001; Begum 
et al. 2006; Kim et al. 2007; Dutta et al. 2008; Roy et al. 
2008; Dutta et al. 2009a, b; Block et al. 2010; Bournaud et 
al. 2010; Dutta et al. 2010; Dutta 2011; Combes et al. 2012; 
Zhang et al. 2012). Note, however, that the uncertainties are 
larger than in the H2 case. For example, high-resolution sim- 
ulations of supersonic turbulence suggest a Burgers scaling 
for both a and E, i.e. a ~ | and b ~ i (e.g., Fleck 1996; 
Kowal & Lazarian 2007; Kowal et al. 2007; Kritsuk et al. 
2007; Schmidt et al. 2008; Price & Federrath 2010). Other 
recent simulation surveys suggest that the scaling exponent 
a is significantly affected by turbulence forcing (Federrath 
et al. 2009, 2010) and self-gravity (Collins et al. 2012). 

In spite of such a burst of interest in ISM turbulence, 
and in spite of the dynamical importance of cold interstel- 
lar gas, there have been very few theoretical works aimed 
at evaluating the effect of turbulence on disc instability. In 
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fact, traditional stability analyses do not take into account 
the scale-dependence of a (or E), but identify a with the typ- 
ical ID velocity dispersion observed at galactic scales. The 
first theoretical work devoted to the gravitational instability 
of turbulent gas discs was made by Elmegreen (1996), who 
assumed Larson-type scaling relations [see Eq. (1)] and in- 
vestigated the case o = — 1 and b — ~. He found that the disc 
is always stable at large scales and unstable at small scales. 
Romeo et al. (2010) also assumed Larson-type scaling rela- 
tions, but explored the whole range of values for a and b. 
They showed that turbulence has an important effect on the 
gravitational instability of the disc: it excites a rich variety 
of stability regimes, several of which have no classical coun- 
terpart. See in particular the 'stability map of turbulence' 
(fig. 1 of Romeo et al. 2010), which illustrates such stability 
regimes and populates them with observations, simulations 
and models of interstellar turbulence. 

In the gravitational instability of galactic discs, there is 
an important interplay between stars and cold interstellar 
gas (e.g., Agertz et al. 2009; Elmegreen 2011; Forbes et al. 
2011; Cacciato et al. 2012). The gravitational coupling be- 
tween these two components does not alter the form of the 
local stability criterion, Q e ff > 1, but makes the effective Q 
parameter different from both the stellar and the gaseous 
Toomre (1964) parameters (Bertin & Romeo 1988; Romeo 
1992, 1994; Elmegreen 1995; Jog 1996; Rafikov 2001; Shen 
& Lou 2003; Elmegreen 2011; Romeo & Wiegert 2011). The 
gravitational coupling between stars and gas also changes 
the least stable wavelength (Jog 1996), among other diag- 
nostics. 

What is the effect of ISM turbulence in this more re- 
alistic context? The first published attempt to answer this 
question was made by Shadmehri & Khajenabi (2012). They 
considered two-component discs of stars and turbulent gas, 
chose a and b so as to sample five of the seven stability 
regimes found by Romeo et al. (2010), and studied the dis- 
persion relation numerically. Their study suggests that tur- 
bulence has a significant effect on disc instability even when 
stars are taken into account. The goal of our paper is to an- 
swer the question above in detail, extending previous work 
along two directions: 

(i) We perform a rigorous stability analysis of two-com- 
ponent turbulent discs, motivated by observations of ISM 
turbulence in nearby galaxies. In particular, we consider two 
complementary cases: Hi plus H2, and gas plus stars. In the 
first case, we examine the dispersion relation analytically, 
and illustrate how the gravitational coupling between H I 
and H2 modifies the main stability regimes of gas turbu- 
lence, which were originally derived neglecting such a cou- 
pling (Romeo et al. 2010). In the second case, we show that 
there are four stability regimes of galactic interest, similar to 
those analysed above, but in only one of them do stars play 
a non-negligible role. We then focus on such a regime, and 
illustrate how gas turbulence affects the onset of gravita- 
tional instability in the disc, i.e. the local stability threshold 
and the corresponding characteristic wavelength. 

(ii) We apply this analysis to a large sample of star-form- 
ing spirals from The Hi Nearby Galaxy Survey (THINGS), 
previously analysed by Leroy et al. (2008) and Romeo & 
Wiegert (2011), and illustrate how ISM turbulence affects 
a full set of stability diagnostics: the condition for star-gas 



decoupling, the effective Q parameter, and the least stable 
wavelength. 

The rest of the paper is organized as follows. The 
(in)stability of two-component turbulent discs is analysed 
in Sect. 2, our application to THINGS spirals is shown in 
Sect. 3, the relation between our results and those of Shad- 
mehri & Khajenabi (2012) is discussed in Sect. 4, and the 
conclusions are drawn in Sect. 5. 



2 (IN)STABILITY OF TWO-COMPONENT 
TURBULENT DISCS 

2.1 Summary of the one-component case 

Here we summarize some of the results found by Romeo et 
al. (2010), which are fundamental to a proper understanding 
of Sects 2.2-2.4. 

The dispersion relation of a turbulent and realistically 
thick gas disc is 

u? = k 2 -2TiGT,(k)k + a 2 (k)k 2 , (2) 

where lo and k are the frequency and the wavenumber of the 
perturbation, and k is the epicyclic frequency. E(fc) and cr(k) 
are the mass column density and the ID velocity dispersion 
measured over a region of size £ — 1/fc, as inferred from 
observations (see, e.g., Elmegreen & Scalo 2004; McKee & 
Ostriker 2007; Romeo et al. 2010): 



£(*) 



* (£)"""• **>-«*(£) 



(3) 



If the disc has volume density p and scale height h, then 
E « 2p£ for £<h and E « 2ph for £>h. The range £<h 
corresponds to the case of 3D turbulence (GMCs and H I 
at small scales) , whereas the range £ > h corresponds to the 
case of 2D turbulence (H I at large scales) . The quantity £0 = 
1/fco introduced in Eq. (3) is the fiducial scale at which E and 
a are observed. This is also the scale at which the Toomre 
parameter Q and other stability quantities are measured, so 
that Qo — ko"o/7tGEo. 

The scaling exponents a and b have an important effect 
on the shape of the dispersion relation [Eq. (2)], and hence 
on the condition for local gravitational instability (co 2 < 0). 
As a and b vary, turbulence drives the disc across seven 
stability regimes, three of which are densely populated by 
observations, simulations and models of galactic turbulence 
(see fig. 1 of Romeo et al. 2010): 

• For b < i (1 + a) and — 2 < a < 1 (hereafter Regime A), 
the stability of the disc is controlled by Qo'- the disc is stable 
at all scales if and only if Qo > Q , where Q depends on 
a, b and £q. This is the domain of Hi turbulence. Both Hi 
observations and high-resolution simulations of supersonic 
turbulence are consistent with the scaling a — b. In such a 
case, the local stability criterion degenerates into Qo > 1, 
as if the disc were non-turbulent and infinitesimally thin. 

• For b > i (1 + a) and —2 < a < 1 (hereafter Regime 
C), the stability of the disc is no longer controlled by Qo: 
the disc is always unstable at small scales (i.e. as k — > 00) 
and stable at large scales (i.e. as k — > 0). 

• For b = i (1 + a) and —2 < a < 1 (hereafter Regime 
B), the disc is in a phase of transition between stability a la 
Toomre (Regime A) and instability at small scales (Regime 
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Figure 1. The main stability regimes of one-component turbulent 
discs (see Sect. 2.1). Also shown are the two-component cases 
illustrated in Fig. 2: Hi in Regime A (a\ = b\ = i) plus H2 in 
Regime A (a 2 = 0.1, b 2 = 0.3), Regime B (a 2 = 0.0, b 2 = 0.5), or 
Regime C (a 2 = -0.1, 62 = 0.7). 



C). This is the domain of H2 turbulence. Note, however, that 
even small deviations from Larson's scaling laws can drive 
the disc into Regime A or Regime C, and thus have a strong 
impact on its gravitational instability. 

Since Regimes A-C are fundamental to a proper un- 
derstanding of Sects 2.2-2.4, we show them in Fig. 1. Note, 
however, that this simple figure is not meant to be a sub- 
stitute for fig. 1 of Romeo et al. (2010), which illustrates all 
seven stability regimes and their relation to the phenomenol- 
ogy of ISM turbulence. 



2.2 Dispersion relation and general properties 

Until now we have considered Hi and H2 separately. How 
does the stability scenario change when H I and H2 are con- 
sidered together? And how does it change when both gas and 
stars are taken into account? We will answer these questions 
here and in Sects 2.3 and 2.4. 

When Hi and H2 are considered together, their gravi- 
tational coupling changes how the disc responds to pertur- 
bations. The dispersion relation can be expressed in a form 
that is particularly useful for discussing the stability prop- 
erties of the disc: 



(u, 2 - M\) (lj 2 - Ml) = (V! - Ml) (Vl - M 2 2 ) 



(4) 



where 



Ml 

,2 



2TvGZi(k)k + a 2 (k)k 2 



(5) 
Vt = K 2 +al(k)k 2 , (6) 

and i — 1, 2p Note that to 2 — M 2 (k) is the one-component 
dispersion relation for potential-density waves [cf. Eq. (2)], 



while V 2 (k) describes sound waves modified by rotation (and 
turbulence). Since M 2 (k) — V 2 (k) represents the self-gravity 
of component i, the right-hand side of Eq. (4) measures the 
strength of gravitational coupling between the two compo- 
nents. 

Eqs (4)-(6) are also applicable to two-component discs 
of gas and stars, even though the stellar component is col- 
lisionless and non-turbulent. This is because stars can be 
accurately treated as a fluid when analysing the stability 
of galactic discs (Bertin & Romeo 1988; Rafikov 2001), and 
because the equations above are valid whether each fluid is 
turbulent or not. Remember, in fact, that the phenomenol- 
ogy of turbulence is encapsulated in Si (k) and <Ji (k) without 
altering the form of those equations. When the disc is made 
of gas (g) and stars (•), S g (fc) and a s (k) are given by Eq. (3), 
while the stellar quantities are not. E*(fc) is the reduced sur- 
face density, E*(fc) = E*/(l + feft*), where the fc-dependent 
factor results from the finite scale height of the stellar layer 
(Vandervoort 1970; Romeo 1992, 1994; Elmegreen 2011). In 
contrast, 0* is the radial velocity dispersion and does not 
depend on fc, since the pressure term in the dispersion rela- 
tion is unaffected by disc thickness (see again Vandervoort 
1970). The gaseous and stellar Toomre parameters are then 
defined as Q g o = fizr g o/7rGE g o and Q+ = k<7*/7tGE*. 

As Eq. (4) is quadratic in lu 2 , it can be solved with 
elementary methods. The discriminant is positive, so there 
are two real roots: 



2 



1[(m 2 + mI)±Va] 



(7) 
(8) 



A = {Ml -Ml) +4(P 2 -X 2 ) (P 2 -M 2 2 ). 

This means that the dispersion relation has two branches 
that do not cross, W+(fc) 7^ uj-(k), except possibly as k — > 
or k — > 00. The functions uj±(k) satisfy two basic properties, 
which constrain the gravitational instability of the disc and 
generalize the stability constraints found in the classical two- 
component case (Jog & Solomon 1984a; Bertin & Romeo 
1988). Such properties are stated and proved below, and 
can easily be visualized with the help of Fig. 2. The cases 
illustrated represent a disc made of marginally stable H I (in 
Regime A) and unstable H2 (in Regimes A-C). 

• Property I: wi(fc) lies below both M\{k) and M 2 {k), 
i.e. a two-component self-gravitating disc is more unstable 
(or less stable) than each component, whether this is turbu- 
lent or not. This can be proved by noting that A is larger 
than (Ml - M 2 ) 2 , so that y/A > \Ml - M%\. In turn, this 
implies that ujI_ < A/f^in, where M m { n is the smallest M 2 
for a given k. 

• Property II: to + (k) is bounded by P\(k) and P 2 (k), 
i.e. this branch is always stable and represents sound waves 
modified by rotation (and turbulence). To prove this, note 
that the inequality \/A > \M 2 — M 2 \ also implies that 
<J\ > A/fmax, where Mm^ is the largest Ml for a given 
k. Note also that lu 2 . cannot be smaller than P^in or larger 
than "Pmaxi otherwise Eq. (4) would not hold. Therefore it 



The dispersion relation of an Af-component turbulent disc is 
y~] _i (M. — P 2 )/(ui 2 — V 2 ) = 1, as can easily be inferred from 
eq. (22) of Rafikov (2001). This equation cannot be expressed in 
a form similar to Eq. (4), and will not be used in the rest of the 
paper. 



2.3 Hi plus H 2 

In Sect. 2.1, we have summarized the main stability regimes 
of one-component turbulent discs. Let us now extend the 
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Figure 2. The two branches of the dispersion relation, iJ\{k) and cj 2 (fc), vs. their one-component counterparts, V 2 {k) and M?(k), 
in stability regimes of galactic interest. These quantities are measured in units of re 2 , the square of the epicyclic frequency, while k is 
measured in units of fcxi,o = re 2 /27rGEio, the Toomre wavenumber of component i = 1 at scale £ = £q. The cases illustrated represent 
a disc made of marginally stable Hi (in Regime A) and unstable H2 (in Regimes A-C). The scaling exponents are specified in Fig. 1. 
The other independent quantities are as follows: fco = 8 fori, (J; Qlfl = 1j £"2,o/Ei,o = ^ anc ^ CT 2,o/o"i,0 = h- These relations imply that 
Q2,o = 2 an< ^ ^0 = 2 ^J2,0> where fcj2,o = 27rGS2,o/o'| is the Jeans wavenumber of component i = 2 (H2) at scale £ = £q. 



discussion to two-component discs of H I and H2 , analysing 
three cases of galactic interest (see again Fig. 2). 



2.3.1 H 



2 in 



nime A 



The response of each component is driven by pressure at 
small scales and by rotation at large scales, while self-gravity 
acts more strongly at intermediate scales (see sect. 2.7 of 
Romeo et al. 2010). This means that the gravitational cou- 
pling between the two components is negligible as k — > and 
k — > 00, and so is the right-hand side of Eq. (4). Therefore 
the two branches of the dispersion relation behave asymp- 
totically as .Mi(fc) and M\(k), i.e. they converge to k 2 as 
k — > and diverge positively as k — > 00. Since the poten- 
tially unstable branch iJLik) lies below Ml(k) (cf. Property 
I) and Mi(k) has a minimum for k > 0, u)_(k) must also 
have a global minimum below n 2 . Thus the disc is stable a 
la Toomre, like each component (see left panel of Fig. 2). 



2.3.2 H 2 in Regime C 

The response of Hi is similar to the previous case, while 
H2 behaves differently (see sect. 2.5 of Romeo et al. 2010). 
The self-gravity term gets dominant for large k and makes 
M'ii2{k) negative. So uj'L^k) is also negative in this limit 
(cf. Property I). For small k, -Mj^W ^ s positive since it is 
dominated by the pressure term (b > 1) and/or the rotation 
term (b < 1). As neither Mii 2 (k) nor M^u^k) is driven by 
self-gravity at large scales, the right-hand side of Eq. (4) is 
negligible as k — > 0. So <JL(k) is positive in this limit, like 



M 2 n (k) and M^ik). The disc is then unstable at small 
scales and stable at large scales, like H2 itself (see right panel 
of Fig. 2). 



2.3.3 H 2 in Regime B 

The behaviour of H2 is intermediate between the previ- 
ous two cases (see sect. 2.3 of Romeo et al. 2010). A sim- 
ilar flow of arguments shows that the disc is in a phase of 
transition between stability a la Toomre and instability at 
small scales, like H2 itself. The middle panel of Fig. 2 il- 
lustrates the phase of small-scale instability, which occurs 
for fco < fcj2,o = 27tGEh2,o/o"h2 ( see Hoffmann 2010 for a 
detailed analysis) . Note how the two components contribute 
to the gravitational instability of the disc, and how their 
coupling widens the range of unstable scales. 



2.4 Gas plus stars 

This case involves three components: Hi, H2 and stars. In 
nearby spiral galaxies, H I and H2 have distinct domains: H I 
dominates the outer regions of the gas disc, while H2 dom- 
inates the inner regions (e.g., Leroy et al. 2008). We can 
then consider Hi and H2 separately. This makes sense here 
because we already know how the gravitational coupling be- 
tween H 1 and H2 modifies the main stability regimes of gas 
turbulence (see Sect. 2.3). What we now want to understand 
is the role that stars play in this stability scenario. Let us 
then distinguish two cases: 
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(i) Stars plus H2. Since the stellar component populates 
Regime A (like Hi) and H2 populates Regimes A-C, this 
case is qualitatively similar to the set of cases analysed in 
Sect. 2.3. So there are three stability regimes: stability a la 
Toomre, instability at small scales, and a phase of stability 
transition. Note that such a variety of regimes is driven by 
H2 turbulence. The stellar component can only modify the 
shape of the dispersion relation; it cannot change the type of 
stability regime. Note also that there is a mismatch between 
two important scales. One is the characteristic scale of stellar 
instabilities, L* = <tJ/7tGE*, which is typically ~ 1 kpc (see, 
e.g., Binney & Tremaine 2008). The other is the largest scale 
at which H2 turbulence has been observed, Lh2 ~ 100 pc 
(e.g., Bolatto et al. 2008). Since L* is one order of magni- 
tude larger than Lh2, the stellar component cannot play a 
significant role in such stability regimes. Therefore this is 
essentially a one-component case, driven and dominated by 
H2. In Sect. 3, we will show that such stability regimes can 
indeed be frequent in nearby star-forming spirals. 

(ii) Stars plus Hi. As both components populate Regime 
A, this is a case of stability a la Toomre: tJL(k) has a global 
minimum, which determines whether the disc is stable for 
all wavenumbers or not (cf. Sect. 2.3.1). In contrast to case 
(i), Hi turbulence reaches scales as large as 1-10 kpc (e.g., 
Kim et al. 2007; Dutta 2011). This makes it possible for 
the stellar component to 'interact' with Hi turbulence and 
contribute significantly to two-fluid instabilities, as in the 
classical case of stars plus non-turbulent gas. 

As discussed above, case (ii) represents the only sta- 
bility regime in which stars play a non-negligible role. We 
then focus on this case, and analyse how gas turbulence af- 
fects the onset of gravitational instability in the disc, i.e. the 
local stability threshold and the corresponding characteris- 
tic wavelength. The effect of disc thickness is well known in 
this context (Romeo 1992, 1994; Elmegreen 2011; Romeo & 
Wiegert 2011). So we do not take that effect into account. 



2.4-1 The stability threshold 

As this is a Toomre-like case, the local stability criterion 
can be expressed in the usual form Q c h > 1, where Q c ff is 
the effective Q parameter. In the classical case of stars plus 
non-turbulent gas, Q e g depends on three parameters: Q+, Q g 
and s — o~ g /o~i,. For analysing Q c a in detail, it is useful to 
factor out the dependence on Q*, Q B a — Q+/Q, and study 
the stability threshold Q as a function of s and q — Q g /Q* 
(Romeo & Wiegert 2011). When gas turbulence is taken into 
account via Eq. (3), Q depends on five parameters: 



10 



so 



OgO 
: 



qo 



a, b and 



Co = £0 fc T * 



(9) 



(10) 



where /ct* = k 2 /2-kGTh, is the stellar Toomre wavenum- 
ber. A general five-parameter study of Q is not more useful 
than a targeted few-parameter analysis. This is because a, b 
and Co are tightly constrained by observations, and because 
their observed values fall within a single stability regime (re- 
member that this is a region of the parameter space where 
the disc has similar stability properties). For these reasons, 
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Figure 3. Contour lines of the stability threshold, Q(so,qo) = 
constant, for discs of stars and turbulent H I vs. the classical case 
of stars plus non-turbulent gas. The solid lines and the shaded 
regions correspond to the median and the la scatter of Co in the 
outer discs of THINGS spirals. 



we analyse Q as a function of so and go, choosing observa- 
tionally motivated values of a, b and Co'- a = b — ~, which 
is the typical scaling of Hi turbulence (see Sect. 1), and 
Co — 1.0 ± 0.5, which are the median and ler scatter of Co 
in the outer discs of THINGS spirals (where Hi dominates; 
see Sect. 3). The range 0.5 < Co < 1.5 is also representative 
of clumpy galaxies at intermediate and high redshiftslj 

Fig. 3 shows a contour map of the stability threshold Q 
for classical and turbulent discs. Consider the classical case 
first, and look at the contour levels Q — 1.1 and Q = 1.4. 
Their slope changes abruptly across the line qo — 1, showing 
that there are two distinct stability regimes. This fact has a 
simple explanation in terms of star-gas decoupling (Bertin 
& Romeo 1988; Romeo & Wiegert 2011). When s < 0.2 and 
go ~ 1, 0J-(k) has two minima: one at small fc, where the re- 
sponse of the stellar component peaks; and the other at large 
k, where gas dominates. For go < 1, the gaseous minimum 
is deeper than the stellar one, and therefore it controls the 
onset of disc instability. Vice versa, for go > 1, it is the stel- 
lar minimum that determines the stability threshold. The 
line go = 1 separates gas- from star-dominated regimes even 
when so > 0.2, but the transition is smooth in this case since 
W_ (k) has a single minimum. In the turbulent case, each con- 



2 Puech (2010) analysed two such galaxy samples at z ai 0.6 and 
zsj2. The median properties of the discs are summarized in his 
table 1 (see also his sect. 3.2). Using those data, we find that the 
stellar Toomre wavenumbcr is fc-r* ~ 0.3 kpc - at z m 0.6 and 
fcx* ~ 0.2 kpc -1 at z ai 2. The spatial resolution is £0 ~ 7 kpc 
and to ~ 5 kpc in the two cases (Puech, private communication). 
This yields Co ~ 2 at z rs 0.6 and Co ~ 1 at z ss 2. Thus, even at 
intermediate and high redshifts, Co is remarkably close to unity 
and lies within the la scatter computed from THINGS. 



© 2012 RAS, MNRAS OOO.rTlfTTI 



6 V. Hoffmann and A. B. Romeo 



10 



OP 1 



0.1 



Turbulent \ 

— — - Classical 



A = 0.1 




0.01 



Figure 4. Contour lines of the characteristic wavelength, 
A(sOiQo) = constant, for discs of stars and turbulent Hi vs. the 
classical case of stars plus non-turbulent gas. The solid lines and 
the shaded regions correspond to the median and the la scatter 
of Co in the outer discs of THINGS spirals. 



tour level is on average shifted down. As Q increases in the 
same direction, this means that turbulence lowers the sta- 
bility threshold, i.e. it tends to stabilize the disc. In Sect. 3, 
we will evaluate the statistical significance of this effect. 



2.4.2 The characteristic wavelength 

The global minimum of w?_(fc) provides another useful sta- 
bility diagnostic: the least stable wavelength A m i n = 27r/fc m i n 
(see Jog 1996 for the classical case). When the disc is 
marginally stable, the value of A m i n is of particular inter- 
est. It is the wavelength at which instability first appears 
as Q* drops below Q. This wavelength can be written as 
A = A At*, where At* = 2n/kT*. The characteristic wave- 
length A depends on the same parameters as Q. So we adopt 
the same approach as before, and analyse A as a function 
of so and qo for observationally motivated values of a, b and 
C a . 

Fig. 4 shows a contour map of the characteristic wave- 
length A for classical and turbulent discs. In the classical 
case, the contour levels A = 0.1 and A = 0.2 are truncated 
above go = 1. This tells us that such short characteristic 
wavelengths occur only when stars and gas are decoupled 
and gas dominates. In fact, star-dominated instabilities ap- 
pear at longer wavelengths: A>0.3 (Bertin & Romeo 1988). 
Note also that the contour A = 0.5 is a separatrix. Levels 
below 0.5 are on the left of this curve (and connected to the 
transition line), while levels above 0.5 are on the right. In 
the turbulent case, each contour level below 0.3 is on av- 
erage shifted to the right, i.e. in the direction of increasing 
A. This means that turbulence shortens the characteristic 
wavelength when stars and gas are decoupled and gas domi- 
nates. An opposite, although weaker, effect is detectable for 



A > 1. Other regimes are also affected, but in a more com- 
plex way. This is especially true for A ~ 0.5, since the sep- 
aratrix of the parameter plane shifts to larger values. Last 
but not least, note how turbulence bends the transition line 
down, favouring star-dominated regimes. In Sect. 3, we will 
analyse these effects in detail. 



3 APPLICATION TO THINGS SPIRALS 

We now consider a sample of twelve nearby star- forming spi- 
rals from THINGS: NGC 628, 3198, 3184, 4736, 3351, 6946, 
3627, 5194, 3521, 2841, 5055, and 7331. For these galaxies, a 
detailed analysis by Leroy et al. (2008) provides high-quality 
measurements of kinematics, as well as stellar and gaseous 
surface densities, at a constant spatial resolution of 800 pc. 

Leroy et al. (2008) also analysed the stability of those 
galaxies, treating the ISM as a single non-turbulent compo- 
nent, gravitationally coupled to stars, with surface density 
Eg = Ehi + Eh2 and velocity dispersion <j g = 11 kms . 
Such a value of <j g fits the H I data well, but is twice as large 
as the typical H2 velocity dispersion observed in nearby spi- 
ral galaxies (Wilson et al. 2011). To represent both Hi and 
H2 well, we choose <r g = 8 kms -1 . This value lies within the 
lcr scatter of ohi (11 ± 3 kms - ; Leroy et al. 2008) and cth2 
(6.1 ± 2.9 kms -1 ; Wilson et al. 2011), and therefore allows 
us to carry out an unbiased stability analysis of THINGS 
spirals. 

The constant spatial resolution of 800 pc used by Leroy 
et al. (2008) makes their data particularly appropriate for 
analysing the effect of Hi turbulence at galactic scales. Hi 
dominates the gas surface density in the outer disc, typically 
for R > 0.43 -R25, where R25 is the optical radius (Leroy et 
al. 2008). We then treat gas as turbulent for R > 0.43 R25, 
and assume Larson-type scaling relations [see Eq. (3)] with 
£0 — 800 pc, E g o = Ego(i?) as tabulated by Leroy et al. 
(2008), and <r g o = 8 kms -1 (see above). Concerning a and 
6, we analyse the case a — b = | in detail, since it represents 
Hi observations fairly well (see Sect. 1). We have also stud- 
ied the case a — b = i as representative of high-resolution 
simulations of supersonic turbulence (see Sect. 1), but here 
we will only mention it when discussing the results of our 
stability analysis. Hereafter we will refer to the model de- 
scribed above as Model 1. 



3.1 The condition for star-gas decoupling 

In Sect. 2.4, we have seen that there is a region in the pa- 
rameter plane where urL(k) has two minima. This is the 
'two-phase region' introduced by Bertin & Romeo (1988) 
and further investigated by Romeo & Wiegert (2011). 

Fig. 5 shows the two-phase region for classical and tur- 
bulent discs. Within this region, stars and gas are dynami- 
cally decoupled and the disc is susceptible to instabilities 
at two different wavelengths, where the responses of the 
two components peak. In the stellar phase the disc is more 
susceptible to long-wavelength instabilities, whereas in the 
gaseous phase it is dominated by short-wavelength instabil- 
ities. Along the transition line between the phases, neither 
component dominates and instabilities occur both at short 
and at long wavelengths. Outside the two-phase region, the 
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two components are strongly coupled and instabilities occur 
at intermediate wavelengths. 

We populate the parameter plane with measurements 
taken from the sample of spiral galaxies, and colour-code 
them by radius. We draw the turbulent two-phase region 
corresponding to the median and la scatter of Co for R > 
0.43 i?25- Note the following points: 

(i) The two-phase region of a classical disc is symmetric 
about q — 1. This symmetry is broken for a turbulent disc 
because gas (dominant for q < 1) follows turbulent scaling, 
but stars (dominant for q > 1) do not. 

(ii) The turbulent two-phase region is larger than the clas- 
sical one. This follows from the fact that turbulence pushes 
the minima of tJL{k) further apart, and the maximum be- 
tween them further up, so as to favour star-gas decoupling. 

(iii) The transition line appears unaffected by the scatter 
of Co- This is because the shape of the two-phase region 
depends on s and q, and q is not affected by turbulence 
(q - qo) if O = b. 

(iv) Turbulence increases the size of the stellar phase 
more than that of the gaseous phase. Recall that the bound- 
ary of the two-phase region is marked by the disappearance 
of the non-dominant peak, i.e. the gas peak in the stellar 
phase and vice versa. Since turbulence affects the gaseous 
peak more than the stellar peak, the size of the stellar 
phase is affected more than that of the gaseous phase. For 
R > 0.43 -R25, this causes a significant number of measure- 
ments to populate the stellar phase. 

(v) For R < 0.43 R25, we find that f 2 = 61% of all points 
populate the two-phase region, two-thirds of them in the 
gaseous phase. In such cases, the onset of gravitational insta- 
bility is controlled by H2. Turbulence is expected to play an 



important role in this process at scales smaller than about 
100 pc (see Sect. 2.4). For R > 0.43 R25, only 4% of all 
points populate this region for a classical disc. This fraction 
increases to 22% for a turbulent disc with a = b — |, and 



to 52% for a = fc= i. 



3.2 The effective Q parameter 

Fig. 6 shows radial profiles of the effective Q parameter, 
Qeff = Qefi(R), for our sample of galaxies. In the left panel, 
we neglect gas turbulence. On the right, we consider turbu- 
lent Hi (a — b = i) for R > 0.43 i?25- Values of Q e ff smaller 
than unity mean gravitational instability. We indicate the 
median and la scatter of Q c s for radii smaller and larger 
than R = 0.43 i?25- We also colour-code the component that 
contributes more to disc instability according to the classi- 
cal condition: gas for Q g o < Q*, and stars for Q* < Q g o 
(Romeo & Wiegert 2011). 

For R < 0.43 -R25, Q B s spans a wide range of val- 
ues, with 13% of points in the unstable regime. Here 56% 
of points are gas-dominated and tend to be less stable 
than the star-dominated points (the median value of Q c s 
is Qeff.g ~ 1.3 and Q e ff,* ~ 2.3 in the two cases). For 
R > 0.43 R25, the range spanned by Q c ff is tighter and only 
4% of measurements are in the unstable regime. Here the 
majority (61%) of points are star-dominated, and there is 
no clear difference in Q c ff between star- and gas-dominated 
points (Qoff, g w 1.3 and Q e ff,* « 1.7). 

Introducing turbulent scaling for R > 0.43 R25 only has 
a small effect on the measurements. For a — b — |, the me- 
dian of Q e ff increases by 3% and the la scatter by 15%. For 
a = b = i the median increases by 6% and the la scatter 
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Figure 6. Radial profile of the effective Q parameter, Q c g(R), for THINGS spirals (Model 1): neglecting gas turbulence (left), and 
taking into account Hi turbulence for R > 0.43 R25 (right). The disc is locally unstable for Q ef j < 1 (light grey shading). For each 
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Thick black lines and dark grey shading indicate the median and 1<t scatter of Q c g in the two radial ranges. 



by 26%. This suggests that turbulence tends to stabilize the 
disc (the median increases), although the magnitude of this 
effect is small and depends on the non-turbulent value of 
Qeff (the scatter increases). 

The stabilizing effect of turbulence seems at odds with 
results from Romeo et al. (2010), who found that the sta- 
bility of gaseous discs is unaffected by turbulence if a — b. 
The difference lies, of course, in the gravitational coupling of 
stars and gas. Consider the approximation for the effective 
Q parameter introduced by Romeo & Wiegert (2011): 



(11) 



r W 1 


if Q* > 


Qcs \ i w 


ifQg> 


w- 2 ^\. 





(12) 



We see that, even if Q g = Q g o, the scaling a g = a g o{£/£o) b 
affects the weight factor W(cT ic ,a s ). The strength of this ef- 
fect is determined by the power-law slope b. Therefore the 
effective Q parameter of turbulent discs always differs from 
the classical case. 



3.3 The least stable wavelength 

Fig. 7 shows radial profiles of the least stable wavelength, 
A m in = Amin(-R), for our sample. On the left we neglect gas 
turbulence, whereas on the right we consider turbulent H I 
for R > 0.43 7?25- Colour-coding indicates the component 
that dominates gravitational instability. As before, the me- 



dian and lcr scatter are indicated separately for small and 
large radii. 

For R < 0.43 R25, there is a clear gap between gas- 
and star-dominated points (the median value of Amin is 
A m in,g ~ 0.7 kpc and A m i n ,* ~ 8.2 kpc in the two cases). 
So the gas-dominated points are characterized by much 
smaller values of A m i n - The discrepancy is less significant 
for R > 0.43 R25, apart from a few measurements close to 
R — 0.43 R25 (A m in,g ~ 3.9 kpc and A m i n ,* « 6.3 kpc). 

Introducing a turbulent gas component for R > 
0.43 R25 causes a significant increase in A m i n . For a = b — i 
the median of A m i n increases by 28% and the la scatter by 
34%. For a = b = |, the median increases by 41% and the 
increase in lcr scatter is again 34%. This suggests a ten- 
dency of turbulence to boost the least stable wavelength. As 
for Qoff , the magnitude of this effect depends on the non- 
turbulent value of A m i n . There is a small number of gas- 
dominated measurements for which the least stable wave- 
length decreases, but these have large uncertainties. 

Why does turbulence affect A c ff more than Q c g? The 
answer is twofold. First, for a purely gaseous disc A m i n in- 
creases markedly with Q c g (Romeo et al. 2010), so that any 
change in Q c g will be amplified in A m i n . Second, as stars 
are taken into account, gas-dominated points can enter the 
star-dominated regime, where A m i n is much larger (see Sect. 
2.4). Both effects depend on the power-law slopes a and b. 
They sum up and drive A m i n to significantly larger values. 



3.4 Robustness of the results 

Modelling the gas disc as a single component with an inter- 
mediate value of a s is not the best that can be done. Here 
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we will no longer follow this traditional approach. We will 
model the gas disc as made of two components, each with 
the more representative value of <j g . A simple way to do it 
is to treat the inner part of the disc as H2 dominated and 
the outer part as Hi dominated. We then set <r g = 6 kms -1 
for R < 0.43 R25 and <r g = 11 kms" 1 for R > 0.43 R25 (cf. 
introductory part of Sect. 3). 

Besides <r g , there is another quantity that deserves 
particular attention: the stellar radial velocity dispersion, 
which we now denote with or*. Leroy et al. (2008) in- 
ferred o\r* from the vertical velocity dispersion, a zir , assum- 
ing that (<7 z /ctr)* = 0.6. In turn, a z ± was inferred from 
the stellar exponential scale height, H+, using the relation 
H* = ct^/2ttG£*. Gerssen & Shapiro Griffin (2012) showed 
that {(Jz/vr)* decreases markedly from early- to late-type 
spirals. The average Hubble stage of THINGS spirals is 
(T) = 4, which corresponds to galaxy type Sbc (the mean 
and the median of T are equal). The best-fitting model of 
Gerssen & Shapiro Griffin (2012) then yields (oz/&r)± = 0.5 
(see their fig. 4). Concerning H+, the relation used by Leroy 
et al. (2008) is not correct. It is the total surface density in 
the disc that determines the stellar exponential scale height: 
H* = al* /27rGE to t, where E to t = E+ + £ g (Bahcall & 
Casertano 1984; Romeo 1992). In view of these facts, we 
set {a z /(7R)i, = 0.5 and use the correct relation for H+. 

Finally, we implement gas turbulence as in Model 1, 
i.e. only for R > 0.43 -R25, where the disc is Hi dominated. 
This is simply because the constant spatial resolution of 800 
pc used by Leroy et al. (2008) is too coarse to probe the 
range of scales affected by H2 turbulence [see Sect. 2.4, case 
(i)]. Hereafter we will refer to the model described above as 
Model 2. 

Table 1 summarizes the dynamical differences between 



Model 2 and Model 1. On the whole, the stability diagnostics 
are moderately affected by the model. The most sensitive 
diagnostic is A m i n , which differs by a factor of 2-3. Q e g is 
more robust, with a difference well below a factor of 2. In 
Model 2, both A m i n and Q B s are smaller for R < 0.43 R25 
and larger for R > 0.43 i?25- 

Despite these differences, the effect of turbulence is com- 
parable in the two models. For R < 0.43 R25, J2 and /2, g 
are slightly larger in Model 2. So H2 is more decoupled from 
stars and slightly more dominant. For R > 0.43 R25, J2 is 
almost identical in the two models, irrespective of the value 
of a = b. Turbulence increases the median value of Q e ff by 
less than 10% in both models. In contrast, the median value 
of A m in increases by 20-30% in Model 2, i.e. less than in 
Model 1. Summarizing, the effect of Hi turbulence in Model 
2 is only slightly weaker than in Model 1. This points to the 
robustness of our results. 



4 DISCUSSION 

Our results cannot be directly compared with those of Shad- 
mehri & Khajenabi (2012), hereafter SK12. This is partly 
because of the wider scope of our paper, which embraces 
a brand-new application to THINGS spirals, and because 
most of the analysis carried out by SKI 2 cannot be easily 
interpreted. 

SK12 analysed five stability regimes of gas turbulence: 
a > 1 and b < | (1 + a); a = 1 and 6^1; and Regimes A-C. 
The first regime corresponds to a fractal dimension D = a+2 
higher than 3, and is therefore beyond the natural range of 
a (see fig. 1 and sect. 3 of Romeo et al. 2010). In the second 
regime, the volume density is scale-independent (D — 3), 
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Table 1. Stability characteristics of THINGS spirals for Models 1 and 2. 



Model 


Radial Range 


a 


b 


/2 a 


/2,g 


/u c 


h d 


Qeff C 


A m in f [kpc] 


1 


R < 0.43 R 25 








61% 


68% 


13% 


56% 


1.67 ±0.68 


1.64 ±2.17 




R > 0.43 R 25 








4% 


100% 


4% 


39% 


1.50 ±0.46 


4.76 ±1.67 






1/3 


1/3 


22% 


41% 


3% 


39% 


1.55 ±0.53 


6.10 ±2.24 






1/2 


1/2 


52% 


9% 


3% 


39% 


1.59 ±0.58 


6.70 ±2.23 


2 


R< 0.43 .R25 








73% 


76% 


25% 


77% 


1.50 ±0.91 


0.67 ±0.62 




R > 0.43 R25 








4% 


100% 


0.5% 


39% 


1.99 ±0.57 


9.32 ±3.40 






1/3 


1/3 


19% 


48% 


0.5% 


39% 


2.09 ± 0.66 


11.10 ±2.78 






1/2 


1/2 


52% 


5% 


0.5% 


39% 


2.14 ±0.72 


12.00 ±2.43 



a Fraction of data that fall within the two-phase region. 
b Fraction of the data points in a that populate the gaseous phase. 
c Fraction of data such that Q c ff < 1. 
d Fraction of data such that Q g o < Q* . 
e Median and lcr scatter of Q c ff. 
Median and lcr scatter of A m ; n . 



so the medium is incompressible and hence subsonic. Cold 
interstellar gas is instead dominated by compressible struc- 
tures and supersonic motions. Therefore even this regime is 
of marginal interest (see again fig. 1 and sect. 3 of Romeo et 
al. 2010). Regimes B and C are populated by H2 turbulence, 
which manifests itself at scales less than Lh2 ~ 100 pc. In 
turn, Lh2 is one order of magnitude smaller than the char- 
acteristic scale of stellar instabilities. Therefore stars play a 
negligible role in these stability regimes [see Sect. 2.4, case 
(i)]. SK12 reached the opposite conclusion. But this is be- 
cause they assumed Larson-type scaling relations even at 
kpc scales, disregarding the type of turbulence associated 
with such regimes. Regime A is populated by both H2 and 
Hi turbulence. While the H2 case raises the same issue as 
Regimes B and C, the Hi case is conceptually simpler. Hi 
turbulence manifests itself at all scales of galactic interest, 
so stars can play a significant role in this stability regime 
[see Sect. 2.4, case (ii)]. SK12 reached a similar conclusion. 
However, even in this case, their approach is different from 
ours. They chose a, b and Co so as to sample Regime A, and 
studied the dispersion relation numerically. We have instead 
examined the whole regime analytically (see Sect. 2.3.1). We 
have then chosen observationally motivated values oi a, b 
and Co, and analysed the onset of gravitational instability 
in the disc (see in particular Sects 2.4.1 and 2.4.2). 

In conclusion, there is a fundamental difference between 
our analysis and that of SK12. Our analysis takes into ac- 
count the astrophysical relevance of the various stability 
regimes, as well as the tight constraints imposed by obser- 
vations of ISM turbulence in the Milky Way and nearby 
galaxies. These are important aspects of the problem, which 
are missing from their analysis. 



(the steeper the Hi scaling relations, the larger the effect). 
This is in contrast to the typical 3-8% increase predicted for 
the effective Q parameter. The effect of H I turbulence is in a 
sense complementary to the effect of disc thickness. In fact, 
disc thickness increases the effective Q parameter by 20- 
50% (Romeo & Wiegert 2011) but hardly changes the least 
stable wavelength (Romeo 1992, 1994) or the condition for 
star-gas decoupling (Romeo & Wiegert 2011). 

Our analysis of THINGS spirals also suggests that H2 
turbulence has a significant effect on the inner regions of 
galactic discs. For _R<0.4i?25, i.e. where H2 dominates 
over Hi, 60-70% of the data fulfil the condition for star- 
gas decoupling and 70-80% of these points represent gas- 
dominated stability regimes. In such cases, the onset of grav- 
itational instability is controlled by H2. Turbulence is ex- 
pected to play an important role in this process at scales 
smaller than about 100 pc (see Sect. 2.4). If a — and 
6 = i then H2 turbulence drives the disc to a regime of 
transition between instability at small scales and stability a 
la Toomre, as was first pointed out by Romeo et al. (2010) 
in the case of one-component turbulent discs. Since this is a 
regime of transition, even small deviations from the standard 
H2 scaling laws (a = and b = i) can have a strong impact 
on the gravitational instability of the disc. This is true even 
when the mass densities of H I and H2 are comparable, since 
small-scale instabilities are more actively controlled by H2 
(see Sect. 2.3). 



5 CONCLUSIONS 

Our analysis of THINGS spirals shows that Hi turbulence 
has a triple effect on the outer regions of galactic discs: (i) it 
weakens the coupling between gas and stars in the develop- 
ment of disc instabilities, (ii) it makes the disc more prone 
to star-dominated than gas-dominated instabilities, and (iii) 
it typically increases the least stable wavelength by 20-40% 
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